Statistical anisotropy of CMB as a probe of 
conformal rolling scenario 

S. R. Ramazanov'^'^*, G. I. Rubtsov^''^^ 

" Niels Bohr International Academy and Discovery Center, Niels Bohr Institute, 
University of Copenhagen, Blegdamsvej 17, DK-2100, Copenhagen, Denmark 
^ Physics Department, Moscow State University, 
Vorobjevy Gory, 119991, Moscow, Russia 
^ Institute for Nuclear Research of the Russian Academy of Sciences 
Prospect of the 60th Anniversary of October 7a, Moscow, Russia, 117312 

May 4, 2012 

Abstract 

Search for the statistical anisotropy in the CMB data is a powerful tool for con- 
straining models of the early Universe. In this paper we focus on the recently proposed 
cosmological scenario with conformal rolling. We consider two sub-scenarios, one of 
which involves a long intermediate stage between conformal rolling and conventional 
hot epoch. Primordial scalar perturbations generated within these sub-scenarios have 
different direction-dependent power spectra, both characterized by a single parameter 
h"^. We search for the signatures of this anisotropy in the seven-year WMAP data using 
quadratic maximum likelihood method, first applied for similar purposes by Hanson 
and Lewis. We confirm the large quadrupole anisotropy detected in V and W bands, 
which has been argued to originate from systematic effects rather than from cosmology. 
We construct an estimator for the parameter h^. In the case of the sub-scenario with 
the intermediate stage we set an upper limit h"^ < 0.045 at the 95% confidence level. 
The constraint on h^ is much weaker in the case of another sub-scenario, where the 
intermediate stage is absent. 
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1 Introduction and summary 



Recent advances in the observational cosmology make it real to start probing the most in- 
triguing aspects of the Universe. In particular, it is of importance to inquire whether the 
statistical isotropy of the scalar perturbations is exact or only approximate. This issue is 
of special interest because the statistical isotropy is one of the key assumptions of the six- 
parametric ACDM model and is favored by inflation. Thus, the violation of this property in 
the observed CMB would imply a highly non-trivial extension of the now standard cosmolog- 
ical model. An additional motivation to search for the statistical anisotropy is the possible 
presence of various anomalies in the CMB data, such as alignment of low multipoles, axis of 
evil, power assymetries, cold spots and others [1]-[12]- 

In the statistically anisotropic but spatially homogeneous Universe, the power spectrum 
of the primordial scalar perturbations ({k.) depends on the direction of the wave vector k. 
The power spectrum can then be written as follows. 



where k = k/k. The coefficients qlai parametrize the direction-dependent part, which 



dependence on the wavenumber k may be absorbed into one function a{k). Commutativity 
of the classical field ((x.) yields Vi^ik) = k) and hence qlm = for odd L. 

The generic prediction of the infiationary theory is that the power spectrum is isotropic, 
a{k) = 0. However, the statistical anisotropy can be generated in models of infiation in- 
volving vector fields [El [151 UHl HZ] or scalar fields with non-minimal kinetic terms [18] ; for 
reviews see also Refs. [19| and [20] . Somewhat more exotic examples are given by introducing 
Bianchi I geometry [21] or noncommutative field theory [221 [23] • The most common feature of 
these models is the statistical anisotropy of a special quadrupole type (the only non- vanishing 
coefficient in a certain reference frame on the celestial sphere is g2o)- This prediction arises, 
e.g., in the model with the rotational invariance broken by a space-like vector [H], or in the 
hybrid infiation incorporating a vector field coupled to the waterfall scalar [15]. However, 
higher multipoles qlm can also emerge within the infiationary framework, see, e.g., Ref. [TS]. 

In infiationary theory, the observed approximate fiatness of the scalar power spectrum 
is due to the approximate de Sitter symmetry of infiating space-time. It has been sug- 
gested some time ago pi] that, alternatively, the fiatness of the power spectrum may be a 
consequence of conformal symmetry. Concrete models of this sort have been proposed in 
Refs. [23 12S] and further developed in Refs. [2Z]^[32]- It is this class of models that we focus 
on in this paper. 

In models of Refs. [23 [26] and similar ones, it is assumed that the cosmological evolution 
starts from or passes through a conformally invariant state with effectively fiat geometry. 





one expands in spherical harmonics YiMik). Unlike in Ref. [T3|, we assume here that the 
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This state is unstable, and conformal symmetry S'0(4,2) gets broken down to 5'0(4, 1) by 
a time-dependent (rolling) scalar field. During this conformal rolling stage, another field 
of zero conformal weight develops perturbations which automatically have flat power spec- 
trum. These perturbations are reprocessed into adiabatic perturbations by one or another 
mechanism (e.g., of Refs. [331 at some later epoch, after the end of conformal rolling. 
It is then natural that the adiabatic perturbations inherit the properties of the original field 
perturbations, modulo possible additional non-Gaussianity. 

At the conformal rolling stage, the properties of perturbations to both linear and lead- 
ing non-linear orders are uniquely determined by the underlying conformal symmetry [32] . 
modulo the overall amplitude and a single dimensionless parameter which we call /i^ (in the 
model of Ref. [25], the amplitude is also determined by /i^, see Section 2.1). This parameter 
governs the non-Gaussianity and statistical anisotropy. The statistical anisotropy generated 
in the conformal rolling scenario is quite different from the predictions of inflation. In par- 
ticular, the coefficients qlm parametrizing the power spectrum ([T]) are the random variables 
rather than fixed parameters of the model. 

There are two sub-scenarios of the conformal rolling scenario which differ by the behaviour 
of the cosmologically interesting modes after the end of the conformal rolling stage. One 
possibility is that these modes are already superhorizon by the end of conformal rolling and 
remain frozen until the late hot epoch. The second sub-scenario assumes that there is an 
intermediate stage at which the field perturbations evolve in a non-trivial way before crossing 
out the cosmological horizon and getting finally frozen out. In the latter case, the evolution 
of the field perturbations results in the statistical anisotropy of all even multipoles. Notably, 
it does not depend on the magnitude of the momentum k, i.e., a{k) is independent of k pJT]. 
The predictions of the former sub-scenario for the statistical anisotropy are considerably 
different. To the linear order in the parameter h, one obtains the statistical anisotropy of 
the general quadrupole type. It is characterized by the amplitude which decreases with 
the wavenumber k, i.e. a{k) oc k~^ [2H]- This means that the corresponding effects in the 
CMB sky are suppressed at high CMB multipoles The statistical anisotropy of the special 
quadrupole type also appears in the next-to-hnear order. Though the effect is suppressed by 
the additional power of the parameter h, now the amplitude a{k) does not depend on the 
wavenumber k [28]. Hence, the subleading contribution may make stronger imprint on the 
CMB. 

Signatures of the statistical anisotropy in the CMB temperature maps have been searched 
for in Refs. [35]-[37]. Motivated by the model of Ref. [H], Groeneboom and Eriksen p5] 
discovered the evidence for the quadrupole statistical anisotropy in the five-year WMAP data. 
However, it was found to be nearly aligned with the ecliptic poles. Using quadratic maximum 
likelihood estimator, Hanson and Lewis [36J extended the analysis to higher multipoles. They 
also included the relevant prefactor in the covariance neglected in Refs. [H] and [35] . Hanson 
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and Lewis [36] confirmed the result on the large quadrupole q2M lying nearly in the ecliptic 
plane. The strongest indication of the statistical isotropy violation, non-zero at the 9(T 
confidence level, was found in the W band of the five-year WMAP data in Ref. These 
findings have been confirmed by the WMAP team [7j in their analysis of the seven-year 
data. One possible explanation of the anomalous quadrupole is the systematics inherent in 
the WMAP data. As argued in Ref. [38], large observed statistical anisotropy may result 
from beam asymmetries rather than have the cosmological origin. 

The purpose of this paper is to constrain the parameter h? of the conformal rolling 
scenario from the non-observation of the cosmological statistical anisotropy in the seven- 
year WMAP data. We follow the general method proposed by Hirata and Seljak for the 
purpose of studying CMB lensing and known as the quadratic maximum likehhood (QML) 
estimation [13]. As discussed in Ref. [36], the same idea can be applied to the study of the 
statistically anisotropic properties of CMB. In this case one assumes that the coefficients 
qLM are small and expands the log-likelihood of the observed CMB to the second order in 
these parameters. By maximizing the log-likelihood with respect to the coefficients qiMi 
one obtains the estimator. Results derived within the QML approximation are in a good 
agreement with the exact likelihood methods. 

We apply this method to construct the estimator for the parameter . In view of the 
results quoted, the estimated values are expected to be inconsistent with the statistical 
isotropy because of the alleged systematics present in the WMAP data. Assuming that the 
interpretation in terms of systematics is correct, we set the upper limits on the parameter 
h? in the following way. For each value of /i^, we simulate the parameter sets {qihi}) and 
then generate a number of anisotropic maps for each set {^la/}- From the maps generated, 
we estimate the values of the parameter h"^. We require that in 95% cases they should not 
exceed the value estimated from the observed CMB. In this way we constrain the conformal 
rolling sub-scenario with the intermediate stage, 

h'^ < 0.045 

at the 95% confidence level. The constraint is much weaker in the framework of the alterna- 
tive sub-scenario. The reason is that the amplitude of the leading order quadrupole decreases 
as k~^. This translates into the suppression of the statistical anisotropy effects at high CMB 
multipoles. Thus, the data useful for the analysis are effectively limited, statistical errors 
are large and the constraint is h? < 190 at the 95% confidence level. The constraint is im- 
proved significantly, once we take into account the subleading contribution to the statistical 
anisotropy. This contribution is of the special quadrupole type, and has the amphtude a{k), 
which is independent of the wavenumber k. Thus, the number of CMB multipoles useful in 
the analysis is much larger. This somewhat compensates the smallness of the constant /i. 
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and we obtain the stronger constraint, 



A 



at the 95% confidence leveL Here Hq is the present value of the Hubble parameter, which 
plays the role of the ultraviolet cutoff, and A is the infrared cutoff. Without going into 
speculation on the value of the constant A, we point out that this constraint is still very 
weak, in view of the fact that the conformal rolling scenario is self-consistent only at /i^ <^ 1 
anyway. 

We conclude that the statistical anisotropy is the relevant signature of the conformal 
rolling with the intermediate stage. It is of particular interest in view of the upcoming 
Planck data. Hopefully, the latter will be free of the quadrupole anomaly. The other expected 
advantage of the Planck data is the larger range of the CMB multipoles, which translates 
into smaller statistical errors. These two factors are expected to improve the sensitivity of 
the data to the parameter h!^ by more than an order of magnitude. On the other hand, 
statistical anisotropy appears to be a weak signature of the alternative sub-scenario, and the 
Planck data are not expected to improve the situation significantly. Thus, it makes sense to 
focus on the other prediction of this sub-scenario, the non-Gaussianity [291 SO] . At the level 
of bispectrum, the non-Gaussianities in the conformal scenario are not particularly special. 
The shape of the intrinsic bispectrum is dictated [39] by the symmetry breaking pattern 
5*0(4,2) — )• S'0(4, 1) and coincides with the bispectrum of a spectator massless scalar field 
in infiationary theory [301 E] (in fact, the intrinsic bispectrum may vanish for symmetry 



reasons, see Section 2.1). The non-Gaussianity generated at the conversion epoch is not 
specific to the conformal scenario either. So, bispectrum alone cannot discriminate between 
the conformal scenario and, say, infiation equipped with the curvaton mechanism. On the 
other hand, the non-Gaussianity of a rather peculiar form arises in the trispectrum. Existing 
constraints (see, e.g. Ref. [12]) are model-dependent and cannot be directly apphed to our 
model. We leave for the future the analysis of the CMB data aiming at the search for the 
non-Gaussianity characteristic of the conformal scenario. 



This paper is organized as follows. In Section 2.1 we review the properties of the field 
perturbations at the conformal rolling stage. The power spectrum of relevant perturbations 
generated by the end of their evolution posseses directional dependence. We show this 



explicitly in Section 2.2 We consider two sub-scenarios: one in which the cosmologically 



interesting modes are superhorizon by the end of the conformal rolling stage (Section 2.2.1) 



and another, with long intermediate stage (Section 2.2.2). In Section [3] we review the main 
ideas of the QML method and construct model-independent estimators for the coefficients 
Qlm- We also construct an estimator for the parameter h"^. Section |4] contains our main 
results. We implement the estimators to the seven-year WMAP data and constrain the 
parameter h"^ of the two sub-scenarios. 
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2 Conformal rolling scenario 

2.1 Fields and their perturbations at conformal rolling 

The main ingredient of the conformal roUing scenario |25l EE] is the conformal stage pre- 
ceding the conventional hot epoch. As we point out below, the properties of the conformal 
rolling stage are quite general [301 122] , as they are unambiguously determined by conformal 
symmetry. Yet it is instructive to review a simple explicit model [22] illustrating this sce- 
nario. It involves a complex scalar field 0, conformally coupled to gravity, which rolls down 
the negative quartic potential V{(j)) = —h'^\(f)\'^, where /i is a small parameter. At large field 
values, the potential is assumed to change and have a minimum at \(f)\ = fo where conformal 
symmetry is explicitly broken. The field is a spectator at the conformal rolling stage and 
somewhat later, i.e., the evolution of the Universe is dominated by some other matter (see 
Ref. [32] for an alternative version where is the dominant matter component). 

The central object in the model is the phase 9 defined by = |0|e*^. Its perturbations 69 
start off as vacuum fluctuations and eventually freeze out. By the end of conformal rolling, 
perturbations 69 have flat power spectrum. Once the radial field |0| settles down to fo and 
conformal symmetry gets broken, what remains are the perturbations of the phase, which at 
this point are isocurvature perturbations. They get reprocessed into adiabatic perturbations 
at much later epoch. In this way the fiat power spectrum of the adiabatic perturbations is 
obtained. It can be slightly tilted if there is small explicit breaking of conformal symmetry 
at the conformal rolling stage fI7\ . 

Let us discuss conformal rolling in more detail. At this stage, the theory is described by 
the action 



where Sq+m is the action for gravity and some matter that dominates the evolution of the 
Universe, and 



is the action for the complex scalar field. Assuming that the background metric is homo- 
geneous, isotropic and spatially fiat, one introduces the field x = where a is the scale 
factor, and obtains its action in conformal coordinates in the Minkowskian form. 



Here rj^'^ is the Minkowski metric and rj is the conformal time. 

Since the scalar potential is negative, the conformally invariant state x = is unstable, 
and the field rolls down its potential. One assumes that the background Xc is spatially 



S — Sg+m + S(j, , 
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homogeneous and without loss of generahty chooses the background solution real. It is given 
by ^ 

h{r]^ - f]) 

where 77^ is a free parameter, which is interpreted as the end-of-roll time. Let us consider 
perturbations in this background. To the leading order in h, perturbations 5xi = \/2SRex 
and Sx2 = V^Imx decouple from each other. We begin with the radial perturbations Sxi- 
They obey the linearized field equation, in momentum representation, 

(6X1)" +p'5xi-Qh\l6xi = 0. 
The properly normalized solution is 



where Bp and are annihilation and creation operators obeying the canonical commu- 
tational relations; H^^^^ is the Hankel function. At late times the solution approaches the 
asymptotics 

= 4^2p,/2^^l_^yBp + h.c. . (3) 

The behaviour 6x1 ~ {v* ~ is interpreted as a shift of the end-of-roll time 77*, which 
now becomes a random field. Indeed, with perturbations included, the radial field Rex = 
Xc + Sxi/ can be written at late times as follows 

Rex = , r . \ T , (4) 

where 

?7^(x) =1]^ + 5r7^(x) , (5a) 

As it stands, the expression Q involves all powers of the small parameter h, i.e., it appears 
to imply some resummation. However, we will be primarily interested in the first non-trivial 
order in h, so we understand Eq. Q merely as a convenient book-keeping tool: to the 
first non-trivial order in h the right hand side of Eq. Q is equivalent to Xc + Sxi- It is 
natural to assume that the fields are initially in their vacuum state. Then the field 5r7^(x) 
is a time-independent Gaussian field with red power spectrum. Clearly, the overall spatially 
homogeneous shift of the end-of-roll time is irrelevant, since it can be absorbed into the 
redefinition of 77*. What is important is the gradient of 77* (x), 

Vi = -9i77^(x) . (6) 
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It has flat power spectrum, while the higher derivatives have blue spectra. 

Let us turn to the object of primary interest, namely, the perturbations of the phase or, 
equivalently, perturbations 8x2 of the imaginary part of the field. We are interested in the 
leading and subleading orders in the small parameter h, so we take the (real) background in 
the form Q. The linearized field equation for 8x2 reads 

{8x2)" - ^i^^8x2 - T^^X2 = . (7) 

At early times, when k{ri^ — ri) ^ 1, we get back to the Minkowskian massless equation. The 
solution to Eq. ([T]) has the following form, 

/d^k (\ 
-— 7r7^('^X2 Hk,x,r7)Ak + /i.c.) , 



(27r)3/V2A; 

where 8x2 ""(k, x, r^) tends to e^^^-'^^'n ^ _oo and Ak, A\. is another set of annihilation 
and creation operators. Modulo corrections proportional to didjrj^, and f^, the solution with 
this initial condition is 



where q = k + kv and v is given by Eq. (|6|. At late times one has 8x2 ^ ['7*(x) — v]~^- 
As a result, the phase perturbations 86 = 8x2/iV^Xi) freeze out. In this regime, the phase 
perturbations, including corrections of order didjr]^ and u^, have the following form |28j : 

where 7 = (1 — is the standard Lorentz factor; and k"^ denote the momenta in 

the direction of the "velocity" v and in the orthogonal direction, respectively. These are the 
phase perturbations by the end of conformal rolling. As promised, they do not depend on 
the conformal time 77, 

dr,8e{x,r])\r,=^^^^) = . (9) 

To the leading order in h, the parameter rj^ is independent of x, so that v = 0, didj-rj^, = 
and the phase perturbations ^ are Gaussian random field with the fiat power spectrum 
and amplitude of order h. To the subleading order, the expression (|8| involves another 
time- independent Gaussian field 5//*. This leads to both non-Gaussianity and statistical 
anisotropy of 86 [x) and hence C(x). Note that the symmetry 6 —9 ensures that the 
intrinsic bispectrum vanishes in the concrete model we review. 

An important remark is in order. Even though we illustrated the conformal rolling 
mechanism by making use of the concrete model ([2]), the results are characteristic of the 
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entire class of conformal models. As an example, the above formulas are valid |30], modulo 
field redefinition, in the Galilean Genesis model [26] based on conformal Galileon field with 
higher derivative action [H]. In fact, these formulas hold [32], provided that the theory has 
the following general properties at the conformal rolling stage: (i) space-time is effectively 
Minkowskian; (ii) conformal symmetry 5*0(4,2) is spontaneously broken down to 5*0(4,1) 
by a homogeneously rolling scalar field of non-zero conformal weight [xi in the above ex- 
ample); (iii) there is another scalar field of zero conformal weight (the phase field 6 above) 
whose perturbations are in the end converted into the adiabatic perturbations. The only 
qualification is that the latter field need not be compact, so the amplitude of its canoni- 
cally normalized perturbations can be arbitrary, whereas the overall amplitude of S6 in the 
model ([2]) is proportional to the same parameter h that determines the amplitude of the 
perturbations Sr]^ and hence non-linear terms in 59, see Eqs. (5b) and 

To avoid confusion, we note that the model ^ is different from the Galilean Genesis 
model [26] and the pseudo-conformal model of Ref. [32] in that in the latter models, the 
cosmological background itself is driven by the rolling scalar field (xi in our notations). 
Therefore, the perturbation of this field (mode 6x1) is gauge-dependent, while the cuvature 
perturbation ( has blue power spectrum. This, however, does not invalidate the discussion 
in the previous paragraph. Indeed, the time shift (analogous to our 67]^:) can be given a 
gauge- invariant definition [26]; this field is related to, but different from (. Furthermore, at 
early times, when the conformal mechanism operates in models of Refs. [261132] . the energy 
density and pressure of the rolling field are small, so the effects due to gravity are negligible. 
This is particularly clear in the Newtonian gauge [2S1 SS] , where the metric perturbations 
are small at early times, while the perturbations in the rolling field have the form equivalent 
to Eq. (g. In this regime, the curvature perturbation ( is irrelevant, as long as one is 
interested in the perturbations of the field similar to our phase 6, and the formulas of this 
Section remain valid. This has been shown expicitly in the pseudo-conformal model of 
Ref. [22], both in the Newtonian gauge and in the gauge 6x1 = 0. 



2.2 Two sub-scenarios and statistical anisotropy of primordial per- 
turbations 

We continue to use the terminology borrowed from the model ([2]). 

At the conformal rolling stage, the phase perturbations 66 get frozen out because of the 
fast evolving background (|4]). Once the conformal rolling stage ends and the radial field 
settles down to a constant value /o, the phase becomes a scalar field minimally coupled 
to gravity. The behavior of its perturbations at later times depends on whether they are 
super- or subhorizon in the conventional sense by the end of conformal rolling. Accordingly, 
there are two sub-scenarios which give different predictions for the statistical anisotropy (and 
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non-Gaussinaity as well) [291 EOl E] . 



2.2.1 Sub-scenario A 

In this sub-scenario, the phase perturbations are already superhorizon by the end of confor- 
mal rolling. So, they remain frozen out, and there is a simple relation between the adiabatic 
and phase perturbations, ({x.) = const ■ 69 (x.) plus possible non-linear terms, where 69 (x.) 
are the phase perturbations late at the conformal rolling stage. So, the properties of ( can 
be read off from Eq. ([s]), modulo possible non-Gaussianity generated when 69 is converted 
into (. 

The interaction of the phase perturbations with the radial ones at the conformal rolling 
stage leads to non-trivial effects in the spectrum of the primordial perturbations. In partic- 
ular, it gives rise to the statistical anisotropy. Indeed, let us consider the two-point product 
69(x)69(x') and average it over the realizations of the operators and A^.. To the leading 
order, we obtain the flat and isotropic power spectrum. The directional dependence appears 
once we take into account corrections coming from the derivatives of the end-of-roll time 
?7=i,(x) and keep only those modes of 6ri^,{x) which are still superhorizon today (shorter modes 
of 5r/=i,(x) give rise to the non-Gaussianity rather than statistical anisotropy [291 EO])- For so 
long modes of 5r7^,(x), it does not make sense to average over the realizations of the oper- 
ators Bp, at this stage. In this way one obtains the power spectrum of the primordial 
perturbations ((k) [2B] : 

Vcik) = Voik) (1 + gi(k) + g2(k)) . (10) 

The directional dependence is encoded in the functions gi(k) and g2(k), which originate 
from the corrections to the linear and next-to-linear orders in the parameter h, respectively, 

gi(k) = -^kkj (^didjr]^, - ^6ijdkdkr]^^ , (11) 

Q,{k) = -^(kv)2 , (12) 

where k = k/k. First, let us consider the leading order contribution gi(k). We expand it 
in spherical harmonics, 

gi(k) = a{k)Y,<lLMYLM{^) , (13) 

LM 

where 

a{k) = k-^ . (14) 



By comparing ( 11 ) with ( 13 ), one concludes that the anisotropic coefficients qlm are Gaussian 



variables, since they are linearly related to the derivatives of the end-of-roll time ?7*(x), which 
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is the Gaussian field. We keep very long modes of St]^,{x.) with p < Hq, where Hq is the 
present value of the Hubble parameter. At shorter wavelengths the field St]^{x.) gets averaged 



out. Therefore, the expression in parenthesis in (11) should be treated as a constant tensor 



throughout our part of the Universe; retaining its dependence on x would result in effects 
suppressed by Ho/k. For this reason, only the quadrupole of the general type survives in 



Eq. (13). Neither its direction nor precise magnitude can be predicted because of the cosmic 
variance. Yet its variance in the ensemble of Universes like ours is calculable. One makes 



use of Eq. ( 5b ) and finds 



{(l2M<l2M') 
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MM' ■ 



(15) 



For similar reason, the second contribution Q2(k) also represents the quadrupole statistical 



anisotropy, but of the special type. It can be expanded in the same fashion as in (13). As 



compared to the previous case, the amplitude a{k) is independent of the wavenumber k. We 
will see that this fact is crucial from the viewpoint of the CMB observations. The other 
important distinction is that the quantities q2M are not Gaussian now. Therefore, it will be 
convenient to work with the components of the "velocity" v, which are Gaussian variables 
with zero means and variances 

3h\ Ho 



In 



(16) 



57r^ A 

Here the present value of the Hubble parameter and the constant A appear as the ultraviolet 
and infrared cutoffs, respectively. The quantities q2M are then given by 



(l2M 



^2m(v) 



(17) 



where v = v/f is the unit vector in the direction of the "velocity" v. 



Equations (14), (15), (16) and (17) are the starting point of our analysis of the statistical 



anisotropy in the CMB within the sub-scenario A. 



2.2.2 Sub-scenario B 

Let us now turn to the second sub-scenario, in which the phase perturbations are still sub- 
horizon at the end of the conformal rolling stage and exit the cosmological horizon later |31| . 
This case is rather non-trivial. First, we need to make certain assumptions about the inter- 
mediate stage following conformal rolling. Barring fine-tuning, we consider the intermediate 
stage very long as compared to the interesting cosmological scales, namely, 

r = k{r]i - r/J > 1 , (18) 

where r]i is the end-time of the intermediate stage, at which the modes S6 exit the cosmo- 
logical horizon. In order not to modify the flat power spectrum generated by the end of 
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conformal rolling, one assumes that the dynamics of the phase perturbations 66 is nearly 
Minkowskian at the intermediate stage, so that they obey the field equation 



dMe = . 



(19) 



The latter condition appears very restrictive. However, there are at least two cosmological 
scenarios where it is obeyed. One is the bouncing Universe with the super-stiff equation of 
state, p ^ p [ini 113 EH] ; interestingly, cosmological contraction with stiff equation of state is 
inherent also in the pseudo-conformal Universe model of Ref. [32]. Another is Genesis |26j . 
where the Universe stays static for a very long period of time before it starts to expand 
rapidly. 

Under the above assumptions, the solution to Eq. (19) supplemented with the initial 
conditions ([s]) and ^ is given by 



h 



47r3/2 7 ^ 



h.c. 



(20) 



where I is the sum of two incoherent waves coming from the direction k and from the opposite 
direction, 

Here 



/ = 1- (^e'^+ [l - kv(+'^) + r{5,, - kikj)drv^'^^'^] + e'^' (^1 - kv^-*^) 



(21) 



ip+ = V'+(x, k) = krji — 2A;?7j,(x + kr), = —krji , 

upper labels (+k) and (— k) denote quanities calculated at the points x + kr and x — kr, re- 
spectively. Note that Eq. (21) is valid at the horizon exit, so it determines, in fact, the prop- 
erties of the adiabatic perturbations (again modulo possible additional non-Gaussianity). 



Note also that under the assumption ( [18| ), the modes of ^r^* relevant in Eq. (21) are indeed 
longer than the size of the visible Universe. 



Proceeding as in Section 2.2.1, we arrive at the power spectrum (10), but now with 

,(+k) 



g(k) = h 



(-k) 



;22) 



In this case all coefficients qlai are generically different from zero. Their variances in an 



(23) 



ensemble of Universes are calculated by making use of Eq. ( 5b ) and are given by [3T 

3/1^ 1 



{(IlmqI'm') 



TT (L-l)(L + 2) 
Notably, the function a{k) does not depend on the momentum, i.e. 

a{k) = 1 . 



(24) 



Equations (23) and (24) determine the statistical anisotropy of the primordial adiabatic 
perturbations in the sub-scenario B in terms of a single unknown parameter h'^. 
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3 Estimators 



3.1 Model-independent analysis 

Let us first apply the quadratic maximum likelihood (QML) method to construct the es- 
timators for the coefficients Qml- Here we define the latter in a model-independent way, 
assuming only that the dependence on the wavenumber k can be factorized as in Eq. ([T]). 
We closely follow the technique developed in Ref. In Section 3.2 we use the same ideas 
when constructing the estimator for the parameter h^. 

In what follows we use the harmonic representation for the temperature fluctuations 
and their covariances unless the opposite stated. The log-likelihood of the observed CMB 
temperature map is given by 

-i:(e|q) = ^etc-^e + ^indetc , (25) 

where q is the vector of coefficients Qlm] the covariance matrix C incorporates the theoretical 
covariance corresponding to the signal as well as the instrumental noise, C = S + N. The 
theoretical covariance is given by 

Here Qim are the theoretical temperature fluctuations of the CMB sky 6T{n) in the harmonic 
representation, 

Gi^ = j dQ6T{n)Y,:^{n) , 

and -P(^(k) is the power spectrum of the primordial perturbations; Ai{k) are transfer func- 
tions. Under the assumption of the statistical anisotropy, we write the theoretical covariance 
as follows, 

S = + 5S , 

where the leading contribution S* comes from the isotropic signal well fitted by the ACDM 
model; the effects of the statistical isotropy violation are incorporated into 6S. The matrix 
S* is diagonal in the harmonic representation, 

Slm;l'm' = Cl6ll'6mm' ■ (26) 

where Ci is the standard CMB angular spectrum. The matrix 6S is given by 

SSi^,i,^> = i''-'Civ (ILM I dn^Y,l^{k)Yv^>{k)YLM{k) , (27) 
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where 



Cu> =4vry d\nkAi{k)Ai,{k)a{k)V^{k) . 
The integral of three spherical harmonics reads 



dn^Y*^{k)Yi,^,{k)YLA,ik) 



1 \m' y^L y^LM 



(2^ 



(29) 



where C^^, are the Clebsch-Gordan coefficients and 



'(2/ + l)(2/' + 1) 



47r(2L + l) 



Normally, the estimators for the coefficients qlm are determined by equating the deriva- 
tive of the log-likelihood to zero, 



dC 







However, the covariance matrix C is not sparse and direct calculations are too costly. Thus, 
we need an appropriate approximation to work with. At this point we make use of the QML 
approach. Assuming that the statistical anisotropy is weak, we expand the log-likelihood 
derivative to the linear order in q. 



dC _ dC 



9qt(9q 



(30) 



We replace the second derivative of the log-likelihood in this expansion by its expectation 
value [36] . 

/ d^C \ _ _ /dC dC 
\ dqdq^ / \dq dq^ 



(31) 



where F is the Fisher matrix. The first equality in Eq. (31) follows from the normalization 
condition 

exp (£) de = l . 



In what follows we use the derivatives of the log-likelihood calculated under the assumption 
of the statistical isotropy unless the opposite stated, and omit the subscript "0" . The first 
derivative of the log-likelihood is then given by 



dC 
9qt 



l0t fey' |g (cy &--tt( (cy ^) 

2 ^ ^ 9qt ^ ^ 2 V 



where C* is the statistically isotropic covariance incorporating the noise, C* = S* + N. The 
second term in the right hand side of this equation is, in fact, the average of the first term 
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over the realizations of CMB. This follows from the identity TrA = (x"''AC^^x), where A 
is any matrix and x is a vector of Gaussian random variables with the covariance C. Thus, 
one writes 

^=h-(h>, (32) 

where 

h^^e.ge, (33) 

and the quantities are the inverse-variance filtered CMB harmonics calculated in the 
absence of the statistical anisotropy, 

e = (S* + N)"^e. (34) 



By substituting Eqs. (31) and (32) into Eq. ([30j) and equating the result to zero, we obtain 
the QML estimator, 

q = (F)-i(h-(h)). (35) 

In what follows we use the Fisher matrix calculated in the full sky and homogeneous noise 
approximation It has only diagonal elements, which do not depend on M, 

2 



v^ (2/ + l)(2r + l) / 

^LM;L'M' — OLL'OmM'W > 



L I I' \ Cii, 



y cf^'cjr' ' 



(36) 



where C/°* is the sum of the standard CMB spectrum Ci and the noise spectrum A^;. The 
constant w denotes the uncut fraction of the sky. We include this factor to achieve better 
agreement between the approximate Fisher matrix and the exact one defined as the average 
over the ensemble of simulated maps with the real sky coverage and inhomogeneous noise. 

3.2 Estimator for the parameter h'^ 

In the framework of the conformal rolling scenario, the coefficients qlm are random variables 
with zero expectation values. To the linear order in h they are Gaussian and have variances 



given by Eq. (15) or (23). The variances depend on the constant h"^, which is the only 
parameter of the model. This makes it possible to constrain the conformal rolling scenario 
from the non-observation of the cosmological statistical anisotropy. We again use the QML 
method to construct the estimator for the parameter h"^. We do that starting from the 
Gaussian hypothesis about the coefficients qLM- This hypothesis is particularly appropriate 
in the context of the sub-scenario B. The non-Gaussian g2Af's appear in the sub-scenario A 
to the subleading order; we comment on this case in the end of this Section. 
The probability density of the coefficients qlm for a given value of h'^ is 
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Here the matrix Q is the covariance of the anisotropy parameters, Qlm;UM' = {qlmQlim')- 
To obtain the hkehhood of the observed CMB with respect to the parameter h"^, one integrates 
the product of two probabihty densities over the set of the parameters Qlm, 

Wi&\h^) = j Wi&\q)Wiq\h^)dq, (37) 



where W(0, q) = exp[£(0, q)] and C is the log-hkehhood introduced in Eq. (25). Following 
the main idea of the QML estimation, we expand the log-likehhood to the second order in 

q, 

where we again replaced the second derivative by its expectation value over the CMB isotropic 



realizations. Now the integral in Eq. (37) takes a simple Gaussian form and can be straight- 
forwardly evaluated, 

, / expf—jF + Q-')-'^] . (38) 

v/det(FQ + 1) '^V2 5qt dqj 



Maximizing (38) with respect to the parameter h'^, 

d\nW{@\h^) _ 

we obtain the equation for the estimator of h"^, 

X.(,PQH-I,-F§)^|^,FQ.I,-g,FQ.I,-f 



In the full sky and homogeneous noise approximation, the Fisher matrix (36) is diagonal, 

FlM;L'M' = Fl^lL'SmM' ■ 

The matrix Q has the same property, 

QlM;L'M' = Qlh^^LL'^MM' , 

where we introduce the quantities Ql which do not depend on the parameter h"^. Then the 
equation determining the estimator takes the form 

(2L + 1)FIQI ^ ^ (2L + l)F,Q, _ 
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Figure 1: Cl of the qLM reconstruction for the V (left) and W (right) bands of the seven-year 
WMAP data. This analysis assumes a{k) = 1 in Eq. ([T|. The la (dark grey) and 2a (light 
grey) confidence intervals are calculated using MC simulated statistically isotropic maps. 
The analysis is performed with the WMAP temperature analysis mask and Imax = 400. 

where we use the same notation h"^ for the estimator as for the parameter of the model. The 



quantities C| entering Eq. ( 39 ) are given by 



M 



where the coefficients qlm are defined by Eq. (35). 



Note that it follows from Eq. (39) that if the predicted statistical anisotropy is of the 
quadrupole form, i.e., with non-zero g2Af's only, then the parameter h"^ can be estimated 
simply as /i^ = C|, modulo obvious additive and multiplicative constants. This is also clear 
on general grounds. Indeed, the rotational invariance requires that the estimator should be 
some function of C|, i.e. h"^ = /(C|). In the small statistical anisotropy approximation, 
we keep only linear terms in the Taylor expansion of the function f{Cl). This immedi- 
ately implies the quoted relationship between h'^ and C|. Since no assumptions about the 
properties of the random quantities q2M have been used in the latter argument, it holds 



for non-Gaussian g2Af's, which describe the quadrupole of the special type, see Eq. (17). 
The only qualification is that the statistical anisotropy is of the order h"^ In ^ in that case. 
Hence, the corresponding estimator reads h^lv? ^ = C|, up to multiplicative and additive 
constants. 
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4 Implementation and results 



We search for the statistical anisotropy using WMAP seven-year maps [50] . We study 
the V and W band data at 61 and 94 GHz. The first step is to implement inverse-variance 



filtering defined by Eq. (34). We write that equation in the form appropriate for applying 



the conjugate gradient technique, 

[(S*)"V Y"fN-^Y]S*e = Y"fN-^0 . (41) 

Here we use the pixel representation for the noise covariance N and the observed CMB 
temperature 0. The matrix Y relates the harmonic space covariance and the observed map, 

where Bi are the beam transfer functions and i labels pixels. We use the foreground reduced 
V and W seven-year maps provided in HEALPix format [32] with Ngide = 512. For 
the beam transfer function we use the average of VI and V2 functions for V band and the 
average of Wl, W2, W3 and WA for W band. 

We consider the noise of the pixels uncorrelated and having the variance a'^/uobs, where 
(To is 3.137 mK and 6.549 mK for V and W bands, respectively, and Uobs is specific to each 
pixel and tabulated in the maps. To remove foreground contaminated pixels we use the 
WMAP temperature analysis mask which leaves us with w = 78% of the sky. We take the 
noise covariance to be infinite (inverse noise is zero) for masked pixels. The noise model 
is constructed using noise covariance and template maps for removing monopole and dipole 
contributions. 

To evaluate the confidence intervals, inverse filtering should be performed on both data 



and large number of simulated maps. Thus, the system (41) must be well preconditioned. 
Following Ref. [36], we make use of the multigrid preconditioner, first proposed by Smith et. 
al. in Ref. [51]. It is known to be the fastest to date and has a typical cost of ten minutes 
when evaluated to lmr,T = 1000. 



Next, we compute the quantities h^M given by Eq. (33). Using Eqs. (27) and (29), we 
write them as follows, 



fiLM - ^ ^(-l)™i' '■Q*^Qrm'Cii>Gjl,C['J^j,_^, . (42) 

lm;l'm' 

We calculate the Clebsch-Gordan coefficients using the GSL [SH] and Slatec [HE] libraries. 



The summation in (36) and (42) is performed up to Imax = 400. We use the publicly available 
Boltzman code (CAMB) [57] to compute the quantities Cw. 

We have simulated large number of statistically isotropic Monte-Carlo (MC) realizations 
of the field using WMAP noise covariance and beam transfer functions. We store the MC 



18 



V-band 



0.035 
0.03 

0.025 
0.02 

0.015 
0.01 

0.005 




+ 



data + 
isotropic 1-0 
isotropic 2-0 



+ + + + + +■ 



2 4 6 8 10 12 14 



0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 




+ 



W-band 



+ 



— 1 1 1 ; 1— 

data + 
isotropic 1 -o 
isotropic 2-0 



+ + + + 



2 4 6 8 10 12 14 



Figure 2: Parameter h? of the sub-scenario B reconstructed from the WMAP V band (left) 
and W band (right). The la (dark grey) and 2a (hght grey) confidence intervals obtained 
from MC simulations are also shown. 

maps in the same format as the original map, and the analysis procedure explained in this 
Section is applied to both data and MC maps on equal footing. 

Now we can check the consistency of the observed CMB with the hypothesis of the 



statistical isotropy. We begin with the model- independent analysis, as outlined in Section 3.1 



We reconstruct coefficients C^, defined by Eq. (40), from the seven-year WMAP data as well 



as from the MC maps. The results are presented in Fig. [T} They are in a good agreement 
with the results obtained by Hanson and Lewis |36] for the five-year maps. In particular, 
we confirm the result on the large quadrupole for the V and W bands. As discussed in 
Refs. [35]-[37], the preferred quadrupole direction lies very close to the ecliptic poles. Another 
suspicious thing is the frequency dependence of the signal. Namely, it is non-zero in the W 
band at much higher confidence level than in the V band. This indicates a systematic 
effect rather than the cosmological origin. As discussed in Ref. |38], the account of beam 
asymmetries can provide a complete explanation of the anomaly. 

Let us turn to the conformal rolling scenario. First, we consider the version of the model 
with the intermediate stage (sub-scenario B). The statistical anisotropy is determined by 



Eqs. (23) and (24). Having the set of the coefficients C| reconstructed from the observed 
CMB, we solve Eq. (39) and estimate the value of h"^. We perform the analysis for the 



multipole numbers starting from Lmin = 2 and ranging up to Lmax = 2 — 14. The results are 
presented in Fig. [2j To evaluate the statistical errors we use about one hundred MC simulated 
isotropic maps. We see that the isotropic model is ruled out at more than 3a confidence 
level even in the V band. However, the large value of h"^ (e.g., h'^ ~ 0.015 at L^ax = 14) 
is due to the anomalous quadrupole anisotropy, which is argued to have non-cosmological 
origin. 
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Figure 3: Parameter h"^ of the sub-scenario B reconstructed from higher multipoles {Lmin = 
4). Results are plotted for the WMAP V band (left) and W band (right). Shown are the 
Icr (dark grey) and 2a (light grey) confidence intervals obtained from MC simulations. 
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Figure 4: C| of the qiu reconstruction for the V (left) and W (right) bands of the seven- 
year WMAP. The momentum dependence of the statistical anisotropy is a{k) = Hok~^. The 
Icr (dark grey) and 2a (light grey) confidence intervals are calculated using MC simulated 
statistically isotropic maps. 



Aiming at constraining the parameter h"^, we simulate a large number of anisotropic maps 
for each value of this parameter. We adapt the approach of Ref . [36] , and use the following 
procedure, adequate in the case of small statistical anisotropy. 



We first simulate a seed map 0* with a covariance S* given by Eq. (26 ). Then we generate 
a set of coefficients {qLAi} based on the value of h"^. The map 



0^ 



I 



6S [S'] ^ 



1/2 



0* 
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has covariance S* + 6S, where SS is given by Eq. (27). To the hnear order in the anisotropic 
effects we have 

Finally, we multiply the map by the beam transfer function in the harmonic space, convert 
it to coordinate space and add pixel noise to get statistically anisotropic simulated map 0*^ 
similar to that observed by WMAP. To set an upper limit, we allow /i^ to be so large that for 
95% of simulated anisotropic maps the value of the estimated parameter exceeds the value 
estimated from the observed CMB map. In this way we obtain the upper limit, which reads 

< 0.045 (43) 

at the 95% confidence level. 

In view of the likely non-cosmological origin of the anomalous quadrupole in the statistical 
anisotropy of the WMAP data, one would like to constrain the parameter from the non- 
observation of higher multipoles only. One way to do that would be to follow the same 



procedure as discussed in Section 3.2 but keep the set {q2M} of the quadrupole coefficients 
fixed and taken from the observational data. In practice, things are simpler. Indeed, the 
effects of the statistical anisotropy corresponding to different multipole numbers L, M do 
not interfere with each other, at least in the approximation of small coefficients Qlm- To see 



this, we note that the theoretical reconstruction of the Fisher matrix (36) is diagonal. The 
covariances of the quantities q^M are diagonal as well. As a consequence, it is straightforward 
to neglect the effect of the quadrupole modulation by using the estimator for the parameter 



h'^ as in (39) but with the summation starting from Lmm = 4. The values of h"^ estimated in 
this way are plotted in Fig. [3j We restrict our analysis to Lmax = 14 and obtain that h'^ is 
consistent with zero for the V band. Making use of the statistical uncertainty inferred from 
isotropic MC maps, we obtain the upper limit 

< 0.040 

at the 95% confidence level. Even though omitting the anomalous quadrupole makes the 



situation cleaner (at least in the V band), this constraint is similar to Eq. (43). The reason 

2L + 1 



is twofold. First, according to Eq. (23), the predicted statistical anisotropy spectrum Cf 
decreases with L as 



^^°^(L-l)(L + 2)- 
Second, the errors grow with the multipole number roughly as L, see Fig. [T] 

With the Planck data available, we expect substantial improvement of the constraint (43). 
The reason is twofold. Hopefully, the quadrupole anomaly will be absent in the Planck data. 
Also, the range of the CMB multipoles useful in the analysis will be considerably extended. 
The error bars, whch can be roughly estimated by making use of the inverse Fisher matrix. 
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scale with the number of muhipoles as /^L- This is clear from the Eq. (36). Taking, e.g., 
Imax = 1200, one would be able to reduce the error bars by about an order of magnitude. 
Hence, the non-observation of the statistical anisotropy will give the constraint as strong as 

< 0.001. We conclude that the statistical anisotropy is a promising signature from the 
viewpoint of the CMB observations in the case of the sub-scenario B. 

Finally, we consider the sub-scenario A. To the linear order in constant h, the statistical 
anisotropy is of the general quadrupole type with decreasing amplitude, a{k) oc k~^. This 
fact is crucial for the search for the statistical anisotropy in the CMB sky. Indeed, the con- 
tribution to the signal 6S is additionally suppressed by the CMB multipole number /. This 



suppression is due to the fact that the integral in Eq. (28) is saturated, roughly speaking, 
at ~ IHq. Effectively, it results in low statistics of the relevant CMB multipoles and large 
statistical errors, which severely restrict the opportunity to observe the (cosmological) statis- 
tical anisotropy of the type predicted. Somewhat loosely we apply the QML estimator to the 
seven-year WMAP data. In Fig. |4] we show the results for of the WMAP reconstructed 
coefficients Qlai, assuming a{k) = Hok~^, but not restricting yet to the quadrupole-only qLM- 
We apply the procedure used in the case of the sub-scenario with the intermediate stage, to 
constrain the sub-scenario A; to this end, the quadrupole point L = 2 in Fig. |4] is relevant 
only. The limit on the parameter h"^ then reads 

< 190 

at the 95% confidence level. Note, however, that for large values of h"^, the QML procedure 
is questionable. This limit can be viewed merely as an indication that the leading order 
contribution to the statistical anisotropy is in fact negligible. The stronger constraint comes 



from the sub leading contribution encoded in (12). The reason is that the amplitude a{k) is 



independent of the wavenumber k in this case. Thus, the suppression at high CMB multipoles 
is absent, and the range of relevant Vs is extended up to Imax = 400. Since the quantities q2M 
are non-Gaussian in this particular case, the constraining procedure is somewhat different. 
First, we generate the components of the "velocity" v starting from a given value of the 
effective constant h'^ln^. Then, using (17), we calculate the coefficients q2M- The quantity 



C| = I '^j^j |g2Mp constructed out of these q2M is compared with the one estimated from 
the seven-year WMAP data. In this way we obtain the constraint, which reads 

hHn^<7 (44) 

at the 95% confidence level. Assuming that the logarithmic enhancement is not particularly 
strong, we conclude that this constraint is still weak. Note also that the statistical anisotropy 
predicted by the sub-scenario A is of the same type as in some infiationary models [H] - [17] . 
Thus, it cannot be used to descriminate our model from other models of the generation of 
primordial perturbations. Fortunately, the sub-scenario A gives rise to the non-Gaussianity 
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in the trispectrum [291 EO] , which is in the sharp contrast with the inflationary predictions. 
We leave for the future search for the corresponding signatures in the CMB sky. 
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